Subnanosecond switching of local spin-exchange coupled to ferromagnets 
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The dynamics of a single spin embedded in a tlie tunnel junction between ferromagnetic contacts 
is strongly affected by the exchange coupling to the tunneling electrons. Using time-dependent 
equation of motion for the spin under influence of the spin-polarized tunneling current, it is shown 
that the magnetic field induced by bias voltage pulses allows for sub-nanosecond switching of the 
local spin and the possibility of spin reversal is illustrated. Furthermore, it is shown that the time- 
evolution of the Larmor frequency sharply peaks around the spin-fiip event, and it is argued that 
this feature can be used as an indicator for the spin-fiip. 
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Detection and manipulation of single spins is an im- 
portant field of science since it pushes the limits of quan- 
tum measurements. Single spins would be an object for 
qubits and would, thus, be crucial for quantum informa- 
tion technology. Potentially, spintronics will replace con- 
ventional electronics devices with spin analogues where 
manipulation, control, and read-out of spins enable func- 
tionality with no or little charge dynamics.^ Of major im- 
portance is to understand how spins can be manipulated 
through electric fields, and questions of the time-scales 
involved. 

Current-induce d magn etic switching have been ad- 
dressed for planai!2l2liIMl and magnetomechanicaP sys- 
tems, as wel l as questions concerning spin-transfer 
^Qj-qy^llIToin] and decoherence.^ Recently, the dynam- 
ics of a local spin coupled to superconducting leads was 
considered and found non-trivial in the sense that a fi- 
nite ch arge current introduces a nutation of the spin 
motion .ESESl Analogous nutations were found for a local 
spin embedded in the tunnel junction between ferromag- 
netic leads biased with a harmonic voltage.^ Concern- 
ing single spin manipulations, such kind of nutation pro- 
vides significant implications for electrically controlling 
the dynamics of local spins. Spin-flip transitions of tun- 
neling electrons generated by spin-orbit coupling is also 
a source for direct manipulations of the local spin.^* Spin 
nutations and spin echo were observed under influence of 
external microwave magnetic field using force detection 
techniques for studies of electron-spin resonance .'i^ 

The purpose of this paper is to consider some implica- 
tions of the nutations of the local spin that are generated 
by a time-dependent bias voltage across the junction.fi^ 
It is argued that that the nutations of the local spin can 
be electrically controlled and its amphtude can be made 
sufficiently large for the spin to undergo spin-flip transi- 
tions, by application of short bias voltage pulses across 
the junction, e.g. pulses of time span less than 1 ns. 
In an arrangement of ferromagnetic leads with unequal 
magnetization and/or in non-coUinear alignment, the lo- 
cal spin tends to line up along the magnetic moment of 
the source lead, and by reversing the bias across the junc- 



tion stimulates reversal of the local spin. In this respect, 
the spin-polarized current can be considered to be gen- 
erating an anisotropy field that stabilize the orientation 
of the local spin, as long as the current fiows across the 
junction. It is further argued that the time-evolution 
of the spin Larmor frequency sharply peaks around the 
event of the spin-flip transition. This feature can thereby 
be used as a signal for detection of the spin-flip event. 

Following the setup given in Ref. H?] a local spin S 
embedded in the tunnel junction between two ferromag- 
netic leads, see Fig. [T] is considered. The leads are di- 
rectly coupled through tunneling and indirectly through 
exchange interaction with the local spin. The Hamilto- 
nian for the system can then be written as 

H = HL + nB + ns + HT- (1) 

The first two terms 7iL(_R) = Y.kaeL{R) ^kacla-^ka de- 
scribe the electrons in the leads, where an electron is cre- 
ated (destroyed) in the left (L) or right (R) lead at the 
energy £ka and spin a =t, I, by c[,^ (cka)- The Hamilto- 
nian for a free spin S in the presence of a magnetic field 
Bq is given by 

ns - -ffMsBo • S, (2) 

where g and /is are the gyromagnetic ratio and Bohr 
magneton, respectively. The two leads are weakly cou- 
pled via the tunneling Hamiltonian 

= X! (4a[^0'^"/3 + TiS ■ aal3]Cqj3 + H.c.) . (3) 
pqajS 




FIG. 1: (Color online) Local spin embedded in the tunnel 
junction between two ferromagnetic leads. 
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Here and henceforth, let electrons in the left (right) lead 
carry the subscript p (q). Here also, aajs is the vec- 
tor of Pauli spin matrices, with spin indices a,f3. The 
rate of the direct tunneling between the leads is denoted 
by To whereas the tunneling with rate Ti is influenced 
by the presence local spin, both of which are assumed 
to be spin-independent. Spin-flip transitions in the di- 
rect tunneling between the leads are neglected in this 
model, and it is furthermore implied that the level of 
the localized electron is far below the Fermi surface of 
the conduction electrons in both electrodes. Therefore, 
the voltage bias should not exceed this energy level dif- 
ference. For convenience the respective amplitudes are 
taken to be momentum independent (although it is not 
required). Typically, from the expansion of the work 
function for tunneling, Ti/Tq ^ J/U, where J and U is 
the spin-spin exchange interaction param eter a nd spin- 
independent tunneling barrier, respectivel}^^^'^ Further, 
an external magnetic field Bo(i) may applied to the sys- 
tem, see Fig. [T] 

It should be noticed that the model should in prin- 
ciple also contain exchange interaction between the lo- 
cal spin and the electrons in the leads, e.g. terms like 
"^Iq^ ■ '^ai3Ck(i, where k & L ov R. Such terms will, how- 
ever, not contribute to the current driven spin dynamics 
considered here, and are therefore omitted. 

Estimates of the signal-to-noise ratio has been pre- 
sented by Balatsky et al^^ and Nussinov et al^^ for 
the present model. It was then found that the spectral 
power density (I^) (^i/^o)^x('^), where x{'^) is a 
Lorentzian associated with power spectrum of the local 
spin. The power spectrum of the shot noise is approxi- 
mately {Ishoti^)) ~ ^0- Using /o ~ 1/Te, where I/t^ is 
the spin-independent scat tering rate. Thus, the signal-to- 
noise ratio becomes {ID / {Ilhoti^)) ^ {Ti/Tofxi^) /re, 
which in the present system is bounded both from above 
and below by a number of order unity.t^Sl Hence, the signal 
may be quite sizable, which enables the further discus- 
sions in this paper. 

When a time-dependent voltage bias is appUed across 
the tunneUng barrier, such that V{t) = Vdc + Vac{t), 
where Vdc and Vac are the dc and ac components, a 
dipole forms around the barrier region through the ac- 
cumulation or depletion of electron charge. This pro- 
cess results in the time dependence of single-particle en- 
ergies as £'p(g) = ep(g) -t- Wi(fl)(t), with the constraint 
WL{t) — WR{t) — eVac{t)- However, the occupation of 
each state in the respective contact remains unchanged 
and is determined by the distribution established before 
the time dependence is turned on. Therefore, the chem- 
ical potentials on the left jiL and on the right lead ijlr 
differs by the dc component of the applied voltage bias, 
Mi ~ A'-R — ^^dc- The tunnel junction is then character- 
ized by two time scales, the Larmor precession frequency 
of the spin lol ~ g/islBj and the characteristic time scale 
of the ac field. 

The dynamics of the local spin was derived in Ref. 
[T71 and was found to satisfy the Landau-Lifshitz-Gilbert 



equatioiP2llSl 



dn . s 



(4) 



under the assumption that the dynamics of the local spin 
is much slower than the electronic processes. This as- 
sumption is motivated by the fact that the energy as- 
sociated with the spin dynamics, hiOL ~ 1 A'-eV, is much 
smaller that the typical electronics energy on the order of 
1 meV.24 In Eq. Q the spin S{t) = Sn{t), where 5 = |SJ, 
whereas n = S/S. The damping factor a{t) oc ST^ /D^^ 
at zero temperature, where 2D is the band widths of the 
leads and it is reasonable to believe this is the case 
also for finite temperatures. Thus, by assuming large D, 
the damping a becomes negligible. 

The effective magnetic field B(t) = Bo(t) -I- B|,j^(i) + 

Bjjj^(t), where the induced magnetic fields can be ex- 
pressed a^ 

5/^sB« (i) = -2ToTi ^ alAfi^P^) " fi^q^)] 

pqa 

xlm /* e'[(?p"-«.-)(*-*')+^(*)-^(*')ldi'z (5) 

J — OO 



gpiB^^it) ^ S f {K(^^{t,t'){ny±-n,y) 

J —OO 

+ [Ki'J{t,t')-Ki'J{t,t')]n.i}dt' (6) 



where 



KgHt,t') = 27rTiReJ2[fiCpa)~mqs 



xe 



i[(«P^-?,s)(t-t')+v(t)-¥'(t')] 



(7) 



and 



KgHt,t')-Kif{t,t')^ 

~27TTiimj2{[mpa) - /(e,^)]e^(«--«'^)(*-*') 



pqa 



-[/(Cp*) -/(e,.)]e^^«''^"«-^^*-*'^}e'[^^*)-'^(*'^l (8) 

Here ip{t) = //^ V{f)dt', for some initial time to, such 
that Lp{t)-ip{t') = Jj! V{t")dt", whereas Cp(,)^ = '^p{q)<y + 

These induced fields are direct responses to the bi- 
asing of the spin-polarized leads. The spin-imbalances 
(^L(i?)T ~ ^L{R)i) in the leads generate a field along the 
z-direction which modifies the Larmor frequency of the 
local spin. More important for the switching, however, 
is whether or not the magnetic moments of the leads 
are equal and/or parallel. Unequal and/or non-parallel 
magnetic moments of the leads generate a magnetic field 
which tends to force the local spin into the magnetic di- 
rection of the source lead, as is shown below. 
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FIG. 2: (Color online) Time-dependent variation of a) 6{t), 
b) jl^ Ki^dt' , and c) the induced Larmor frequency ujY^'^{t), 
for a square bias voltage pulse, c.f. dotted line in a). Here, 
Vac = 1 mV, TiNo = 0.1, PL = -PR = 1/2. 



In order to better understand the motion of the lo- 
cal spin, consider the corresponding classical equation 
of motion, which is consistent with the parametriza- 
tion of the spin on the unit sphere, i.e. Sn — 
S {cos (/) sin 6, sin (f> sin 9, cos 6). With Bq = Bqz time- 
independent, the equations can be written 



+S [ [Ki^J{t,t') - KfXt,t')]dt' cose, 



dt 



dd 



(9) 



^ S j Kg\t,t')dt' sine. 

J —OQ 



The equation for the polar angle e can be analytically 
solved and in the stationary regime it gives 



9{t) — 2 arctan 



tan^e-2-^-^?^= 
2 



,NLaNRs{t-to) 



(10) 

where = ^(^o)- This expression clearly shows the pos- 
sibility to switch the local spin by means of a stationary 
bias voltage. For example, configuring the leads such 
that Nl^Nui — iVij^iVfli > O,'^ and biasing the system 
by Vdc > 0, gives an exponential decay of the argument 
which leads to that the polar angle e{t) (local spin 
becomes 1), as t — > oo, for any given initial polar angle 
eQ. By the same argument e{t) -k (local spin becomes 
I) if the bias Vdc < is applied. 

Using the result for the polar angle within the sta- 
tionary regime, the characteristic time-scale K/tc ~ 
2TTeVTl (Jaa^LaNRa for the polar angle motion of 
the spin is obtained. Parametrizing the spin-polarized 



DOS NLiR)a = No{l+apLiR))/2^wheK -1 < pl(H) < 
1, assuming TiNq ~ 0.1, pl = —PR = —1/2, and V ^ I 
mV gives the characteristic time scale r^, « 5 ps, which 
is sufficiently short to switch the spin from being J, to | 
within 1 ns, see Fig. |2]a). 

In a realistic set-up of the system it would be desirable 
to obtain the current induced switching by application of 
bias pulses of some time span . A sudden onset of a bias 
pulse leads to transient effects in the induced magnetic 
fields B|j"j(t), n = 1,2, which transfer into the motion 
of the spin. Therefore, consider application of the bias 
Vrfc+Kc[B(t-To)-6(t-ri)], where n to+t^. We also 
assume that Bq = since we are interested in the local 
spin dynamics generated by the spin-polarized current. 
As shown above, a stationary bias will eventually line up 
the local spin in the magnetic direction of the source lead 
and since the main interest lies in a switching obtained 
by a pulsed bias, henceforth Vdc = 0. Then, the equation 
of motion for e(t) can be written 



de 

lit 



x[e(i-ro)-e(i-ri)]+ e-(*--i)/- (11) 
cos [eVacit - Ti)]e(i - ri)l sine 



-(t-To)/r 



where the time-scale r relates to the electronic tunneling 
processes. Such processes are of the order of 1 fs, which is 
much faster than the characteristic time-scale Tc and the 
Larmor frequency Hlul for biases between 1 - 100 mV.'^ 
For e.g. a 1 ns bias pulse of 1 mV, the critical time scale 
is 100 fs, which may be achieved within present experi- 
mental state-of-the-art-technology for nanoscale systems. 
Physically the time-scale r implies that the induced mag- 
netic field is a retarded response to the bias across the 
junction. 

The calculated time-dependence of the polar angle is 
plotted in Fig. |2] a) for a square bias voltage pulse, and 
it is readily seen that the angle goes from tt to within 
the time span of the pulse, hence, the pulse is sufficiently 
long to ffip the spin from j to f. The plot in Fig. [2] 
b) shows the induced magnetic ffeld which is acting to 
align the local spin along the magnetic direction of the 
source lead. At the onset (termination) of the bias pulse, 
the amplitude of the induced magnetic field grows (de- 
cays) exponentially, as expected from Eq. (111. At the 



pulse termination, however, there are additional oscilla- 
tions in the induced field, as a reaction to the removed 
bias. These oscillations, are not visible in the motion of 
the polar angle, since they are exponentially suppressed. 



c.f. Eq. (10 1 



The Larmor frequency of the precession of the lo- 
cal spin is affected by the time-dependent variation of e 
that gives rise to a momentary change of its value. By 
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t-to(ns) 



FIG. 3: (Color online) Time-dependent variation of a) 9{t) 
and b) ui^£'^{t), for different initial polar angles 60, under the 
bias V{t) = Fdc + Vac{[Q{t - to) - e{t - ti)] - [e{t - ti) - 
Q{t — t2)]}, with — to = ^2 — ii = 1 US. For clarity, the plots 
in panel b) have been shifted 5 GHz upwards for increasing 
9o. Other parameters are as in Fig. [2] 



definition 



u;l{t)^{g^^sr\B{t)\'={B„+Bl^l{t) 



+S / [Ki'J{t,t') - Kif{t,t')]dt' cos e{t) 

+ (^sj K^/Jit,t')dt'^ sm'eit) (12) 

The time-dependent parts of the induced fields in the 
first term are equal and can be written as (Vdc = 0) 



eVa, 



-{t-r„)/T _ -(t-ri)/T 



sin [eVac{t-Ti)]Q{t-n). 

(13) 

Hence, this field does not affect the Larmor frequency 
until after the bias pulse has terminated, and then it gen- 
erates an oscillatory variation of wz, with an exponential 
decay, see Fig. |2]c). These oscillations are sufficiently 
large to enable read-out of the time instant when the 
pulse is terminated. 
The last term in Eq. 



(12 



containing the field that 
generates the spin-fiip, vanishes as long as 9(t) = 0, ±7r, 
that is, when the spin is either | or J, in the global spin 
quantization axis. This term gives a contribution, how- 
ever, during the short time interval when the spin flips, 
since then < \9{t)\ < tt. Hence, one would expect a 
sharp peak in the time-evolution of LOL(t), which indeed 
is seen in Fig. |2]c). This peak in the time-evolution of the 
Larmor frequency is a fingerprint of an actual occurrence 



of a spin-flip and would be measurable. Using the values 
comparable to those in Fig. [2] results in an amplitude of 
the peak of the order of 50-500 GHz in case of vanishing 
external magnetic fleld, for biases 1-10 mV. 

In the discussion we have neglected influences from e.g. 
anisotropy fields. The absence of such or similar fields 
that would stabilize the spin direction, cause the spin 
orientation to drift randomly in equilibrium. It has been 
shown that the spin-polarized current between two ferro- 
magnetic leads has such effect on the local spin, i.e. the 
effect of the spin-polarized current on the local spin is to 
stabilize its orientation. The plots in Fig. [sja) show 9{t) 
for several different initial polar angles 6*0, given the bias 

vit) = + Vacimt - to) - e{t - 1,)] - [e{t - 1,) - 

0(i - *2)]}, with ti - to = t2 - ti ^ 1 ns. The plots 
illustrate how the spin-polarized current forces the local 
spin to be aligned with the majority spin of the source 
lead, 9(t) — s- for tg < t < ti. Furthermore, by reversing 
the bias voltage, the spin orientation can be reversed at 
will, which is illustrated by the plots in Fig. |3] a) where 
9{t) TT, in the interval ti < t < t2- Panel b) in Fig. 
[3] further illustrates how the Larmor frequency peaks at 
the instant of the spin-flip, both when its polar angle ap- 
proaches zero and tt. The increasing amplitude of the 
peak with increasing difference between the final and ini- 
tial polar angle is expected, because of the increasing 
reorientation path the spin has to traverse. 

The result displayed in Eq. (Ill shows that it would 



be possible to fiip the spin by having only one of the leads 
ferromagnetic and the other non-magnetic. Thus, a set- 
up of an scanning tunneling microscope (STM) with a 
spin-polarized STM tip and non-magnetic substrate sur- 
face would correspond to having, say, pL and pn — O. 
Hence, letting pL = 1/2 and using the same parameters 
as in Fig. [2] one finds a time scale for the spin flip of 
roughly ~ 10 ps, that is, doubled the time scale com- 
pared to the case discussed above. The spin flip effect 
should nonetheless be likewise observable as in the case of 
two ferromagnetic leads in anti-parallel alignment, how- 
ever, on a larger time scale. 

The result in Eq. ( 11 1 also shows that if both leads 
are being ferromagnetic and in parallel alignment, the 
resulting spin flip effect on the local spin decreases since 
the induced fields from the left and right leads tend to 
cancel each other. The polar angle becomes a constant of 
motion whenever the two leads are equally strong ferro- 
magnets in parallel configuration, since then p^ = pr = p 
which yields KaNLaNR^ = (7Vo/2)2[(l +p)(l ~p)- 
(l-p)(l-Hp)] = 0. 

Anisotropy field acting on the local spin has been omit- 
ted in the present paper for simplicity. In a realistic sys- 
tem, however, the local spin is likely subject to some 
type of anisotropy field which stabilize definite spin di- 
rections. Such fields may be small enough to not substan- 
tially increase the time scale of the spin-fiip processes in 
the present situation. Without any type of anisotropic 
field acting on the local spin, its orientation will in equi- 
librium, i.e. V{t) = 0, most likely describe a random drift 
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due to exchange interaction with equiHbrium fluctuations 
in the current between the leads. Such a current will lack 
any definite spin-polarization that the non-equilibrium 
situation can present. 

In summary, it is demonstrated that the dynamics of a 
local spin embedded in the tunnel junction between fer- 
romagnetic leads can be manipulated at will by means 
of a time-dependent bias voltage. Especially, spin-flip 
transitions can be stimulated at a sub-nanosecond time- 
scale by bias voltage pulses, even for a moderate spin- 



polarization of the leads. The time-evolution of the Lar- 
mor frequency of the local spin sharply peaks at the spin- 
flip event, something that determines actual occurrence 
of a spin-flip. These flndings would be useful in the ad- 
vent spintronics and quantum information technology. 
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